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Abstract 

We present the results of Monte Carlo simulations of a model of a single maghemite ferrimagnetic nanoparticle 
with the aim to clarify the role played by the increased anisotropy at the surface and by the shape (spherical 
or elliptical) of the particle on the magnetization processes at low temperatures. The formation of hedgehog-like 
structures for high enough surface anisotropy is responsible for a change in the reversal mechanism of the particles. 
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The progressive size reduction of magnetic par- 
ticles with technological application has brought 
renewed interest on the issue of surface and finite- 
size effects in their magnetic properties. In partic- 
ular, it is an old-known fact [1] that, when entering 
the nanometric range, the reduction of lattice sym- 
metry at the surface of the particle results in an 
enhancement of magnetic surface anisotropy with 
respect to bulk values which is not clearly under- 
stood at present, although a number of computer 
simulations have been recently performed in order 
to clarify this issue [2]. Here we present the re- 
sults of Monte Carlo simulations of a model of a 
ferrimagnetic nanoparticle with the spinel lattice 
structure of 7— Fe203 (maghemite). We consider 
an extension of our previous model in [3] to the 
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Heisenberg case, including the term i/anis/^s — 
- Eti {KciSff + Ks{S, ■ n,)2) to account for 
finite anisotropy. Here surface spins are allowed to 
have radial surface anisotropy Ks distinct from the 
uniaxial anisotropy Kq of those in the core. 

By using the standard Metropolis algorithm for 
continuous spins [5], we have computed the spin 
configurations of ellipsoidal particles attained after 
a cooling from a high temperature phase. We have 
found that, for high enough Ks/Kc ratios, these 
configurations change from an AF state in which 
spins in each sublattice are almost aligned along 
the core easy-axis to a hedgehog-like structure with 
spins pointing along radial direction and AF order 
among sublattices (not shown here due to space 
limitations [6]). These structures are not observed 
when uniaxial surface anisotropy is considered and 
therefore can be attributed to the increased surface 
anisotropy. 



Preprint submitted to Elsevier Science 



2 February 2008 



-200 -100 100 200 
h(K) 

Fig. 1. Hysteresis loops for a spherical particle with di- 
ameter D = 3a, Kc = 1, (a) Ks = 10, (b) Ks = 50. 
The main panels show the total magnetization along the 
z-axis (circles) and the corresponding contribution of the 
core (dashed lines). The upper subpanels show the surface 
(circles) and core (continuous lines) contribution to M„ 
defined in the text. 

Zero temperature hysteresis loops for a spherical 
particle with diameter D = 3a and an ellipsoidal 
particle with the same diameter and major axis 
L = 8a with Kc = 1 and Ks = 10, 50 are com- 
pared in Figs. 1 and 2. In these figures, we also show 
the sum of projections surface and core magne- 
tizations into the respective anisotropy directions 
Mn — ^iLi I'^i ' First, we notice that as parti- 
cles become more elongated, although the coercive 
field and remanence magnetization remain almost 
constant, the hysteresis loops become more elon- 
gated ressembling those observed experimentally. 
Also the reversal of the core (dashed line) becomes 
less uniform, a fact that can be attributed to the 
increasing proportion of spins with anisotropy di- 
rection far from the field direction when increas- 
ing L. As Ks increases from 10 to 50 the coercive 
field he increases and, more importantly, there is 
a change in the magnetization reversal mechanism 
both for spherical and elliptical particles. Whereas 
in the first case the particle core reverses in a quasi- 
uniform form with the spins pointing mostly along 
the z-axis {M„ « 1) during the process except near 
he and with the surface spins following the core re- 
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Fig. 2. Same as in Fig. 1 but for an elliptical particle with 
D = 3a, and long axis L = 8a. 

versal (the departure ). However, at higher Ks [(b) 
panels], surface spins remain close to the local ra- 
dial direction (M„ w 1) during all the reversal pro- 
cess, driving the core spins out of their local easy 
axis and making their reversal non-uniform (0.5 < 
Mn < 1) due to the appearance of the hedgehog- 
like structures mentioned above. 
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